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Abstract. When written in terms of ?9- functions, the classical Frobenius Stickclbcrgcr 
pseudo-addition formula takes a very simple form. Generalizations of this functional identity 
are studied, where the functions involved arc derivatives (including derivatives with respect 
to the modular parameter) of the elliptic trilogarithm function introduced by Bcilinson and 
Levin. A differential identity satisfied by this function is also derived. These generalized 
Frobenius-Stickelberger identities play a fundamental role in the development of elliptic 
solutions of the Witten-Dijkgraaf-Verlinde-Verlinde equations of associativity, with the 
simplest case reducing to the above mentioned differential identity. 
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1 Introduction 



Amongst the many beautiful identities satisfied by the elliptic (and related) functions is the 
Frobenius-Stickelberger relation [12} [22] 

(C(a) + ((b) + C(c)) 2 = p(a) + p(b) + p(c), 



(a + b + c = 0). 
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Writing this identity not in terms of Weierstrass functions but in terms of ^-functions yields the 
equivalent form 



■&[(c) 0j(a) 

01 (a) tfi (6) ' 0i (6) 0i (c) ' 0i (c) 0i (a) 



0i(a) 0^(6) + 0i(&) 0i(c) 



+ 



+ 



1 



0'/(a) + + *jf(c) 



01 (a) Mb) 01 (c) 



1 0f (0) 

2 0^(0) ' 



(2) 



where again a + b + c = 0. With the help of the heat equation the second set of terms may be 
written in terms of derivatives with respect to the modular parameter r. 

The purpose of this paper is to explore certain neo-classical identities satisfied by the elliptic 
trilogarithm introduced by Beilinson and Levin [HE]. This function will be defined in Section[2l 
but for now it is sufficient to note that with it the above identity takes the simplified form 

/( 3 '°)(a)/( 3 '°)(6) 

+ / (3,0) (6)/ (3,0) (c) 

I +/( 3 > >(c)/( 3 >°)(a) 
These new identities take the schematic form 

{quadratic terms} + {linear terms} = {0 constants} , (3) 



{f^\a) + f^(b) + f^(c)} = 



r(2,l) 



(2,1), 



*This paper is a contribution to the Proceedings of the Workshop "Elliptic Integrable Systems, Isomonodromy 
Problems, and Hypergeometric Functions" (July 21-25, 2008, MPIM, Bonn, Germany). The full collection is 
available at |http://www. emis.de/journals/SIGMA/Elliptic- Integrable-Systems.html] 
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where the linear terms contain one more r-derivative than the total number of r-derivatives in 
each part in the quadratic term. Before these identities are discussed (Section HJ a differential 
equation satisfied by this function f(z,r) will be derived in Section [3j Applications of these 
new identities are then discussed in Section We begin by defining the elliptic polylogarithm. 



2 The elliptic poly logarithms 

The classical polylogarithm is defined, for \z\ < 1, by 



oo r 



, n' 

71=1 



and by analytic continuation elsewhere. A first attempt at an elliptic analogue of this function 
might be 

oo 

n r (C,q)= £ Li r (g"C)- 

n=— oo 

However this series diverges, but by using the inversion formula (jlip and ^-function regulariza- 
tion one can arrive at the following definition of the elliptic polylogarithm function [Tl I14j: 

oo oo 

£i7-(C,Q) = £Lv(g n C) + £Li r (^r 1 ) -Xr(C,?), r odd, 

n=0 n=l 

where 

xAc,,) -E (r _^. +1) , aogc) ( - j) (iog^. 

A real- valued version of this function had previously been studied by Zagier [23] . With this the 
function / may be defined. 

Definition 1. The function f(z,r), where z G C, r G H, is defined to be: 
/(.,r) = ^{£i 3 ( e 2 ^, g )-£i3(l, g )}. 
The function /( n ' m ) that appear in the introduction are the derivatives 

f(n,m) J_ 

~ dz n dr m ' 

It immediately follows from the definition that 
d \ 3 1 x 1 



£i 3 (l,g) = — E a (t). (4) 



and 



dry (2vri) 3 ° v 120 



9 V 1 f^i(^r) 

/(z,r) = -— log 



9z y 27ri [ v( T ) 

Thus the elliptic-trilogarithm may be thought of as a classical function (or, at least, a neoclassical 
function) as it may be obtained from classical elliptic functions via nested integration and other 
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standard procedures. It does, however, provide a systematic way to deal with the arbitrary 
functions that would appear this way. 

The following proposition describes the fundamental transformation properties of the func- 
tion: these will be used in subsequent section to prove the various differential and functional 
identities. The precise definitions and normalizations of the various objects used here are given 
in Appendix A. Also, the notation F ~ G will be used if the functions F and G differ by 
a quadratic function in the variables z and r. These quadratic terms may be easily derived, but 
will play no part in the rest of the paper. 

Proposition 1. The function f has the following transformation properties: 



f(z + l,T)~f(z,T); 
/(z,r + l) = /(z,r); 

f(z + r, r) ~ f(z, t) + | ^z 3 + ^z 2 t + ^zt 2 + ^r 3 
f(-z,r) ~ f(z,r). 
The function also has the alternative expansions: 
1 f 1 



(2iri) 



-z logz + z lo gV (r)j +J —^^ (2n + 2)!(2n ) (^) 



(5) 



and 



sin 2 (7rrz) 



(6) 



XTie ,/srsi o/ i/iese imply the following the transformation property: 



If! 

r 3 4 ! ■ 



Proof. The first three relation follow immediately from the definition. The fourth used the 
inversion formula for poly logarithms pip. 

The proof of (J5J) and (jSJ) just involves some careful resumming. Consider the first two terms 
in the definition of /: 



j^Lia^e**-) +X> 3 (<fe- 27r -) = Li 3 (e 2 ™) +2^ 



n=0 



n=l 



s=0 



i-ir 

(2«)! 



E 1 

n,r=l 



raTy2s-3 



(27T2) 



2.s 



From this series © follows immediately. To obtain ([5]) one rearranges the terms. The s = 
term cancels in the final expression and the remaining terms may be re-expressed in terms of 
Eisenstein series (for s > 1) or the Dedekind function (for s = 1). Finally, using the result 



^ ? ^u,(^) = -i[i + «*h(™)] 



one may obtain a series for L13 (e 2mz ) . Putting all these parts together gives the series (J5J). 
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(z + t) 


= f (m (z 




(z + t) 






(z + t) 






(z + t) 


= f m {z 



3 Differential identities 

Theorem 1. The function f satisfies the equation 

/ (3 ! 0) / (l ! 2)_ (/ (2 ! l) ) 2 + l / (0 > 3) = __l_^ (r) 

Proof. We denote the left hand side of the differential equation by A(/)(z,r) and study its 
transformation properties. It follows from Proposition [1] that the third derivatives of / are 
invariant under the transformation z^ z+ \ and that under the transformation z i— ► z + r one 
has: 



f {m (z)-\, 
2f^(z) + f^\z) + \, 
3f^(z)+3f^\z)-f^(z)- 1 -. 

It immediately follows that combination A(/) is a doubly periodic function (even though 
individual components are not). From the Laurent expansions of /, the only term with a pole 
is f( s '°\ which has a simple pole at z = 0, but this cancels with the zero at z = of the term f^' 2 \ 
Thus A(/) is a doubly period function with no poles, and hence must be independent of z, i.e. 
a ^-constant. 

Using the series expansion in Proposition [1] one finds that 

A = -{— E'--E\ 
12 \2vri 2 12 2 

and hence the result follows on using the Ramanujan identity 

0IE2 E 2 — E\ 
q ~dq~ = 12 ' 

An alternative proof is to study the modularity properties of A under r — ► — r^ 1 . One finds 
that A must be a modular form of degree 4, and hence must be a multiple of E4, the constant of 
proportionality being straightforward to calculated. Note, one could easily redefine /, using 
so that A = 0. ■ 

This differential identity contains much information. Using the z-expansion of the function / 
it is equivalent to an infinite family of identities between the Eisenstein series and their deriva- 
tives, though these will all be consequences of the Ramanujan relations and the fact that the 
ring of modular forms is finitely generated by E4 and Eq. Alternatively, using the (/-expansion 
of / it is equivalent to an infinite family of number-theoretic trigonometric identities (number 
theoretic since in the g-series for / the q n coefficient will involve a sum over r\n). 



4 Functional identities and rank 2-root systems 

The basic identity ([2]) may be associated to the ^-root system via the identification a = (a, z), 
b = ((3,z), c = —(a + (3,z), where a and (3 are the positive roots and (•,•) is the standard 
Euclidean inner product. This immediately suggest that there should be variants of this identity 
for other rank-2 root systems. 
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Lemma 1. Let 1Z be the root system for the 2- dimensional Coxeter groups A 2 , B 2 or G 2 , with 
the standard normalization for a, (3 positive simple roots: 



A 2 
G 2 



(a,a) = ((3,(3)=2, (a,/?) = -l, 
(a,a)=2, (P,P) = 1, (a, (3) = -1, 
(a, a) =6, {(3,(3) = 2, (a,/3) = -3. 



Then 



^ (a, /3)/ (3 ' 0) ((z, a), r) • / 3 '°) ((z, (3), r) + £ k a f^ ((z, a), r) = 0, 



where: 



A 2 
B 2 
G 2 



k a = 1 /or a// roots; 



^short 2, 
^short = 10, 



A- 



long 



k 



long 



6. 



The proof is entirely standard and is omitted: it follows the proof of Theorem [TJ one first 
shows that the left hand side is doubly periodic in all a, b, c variables (constrained via a+b+c = 0) 
with no poles in all these variables, and hence must be independent of all three variables, i.e. it 
must be a ^-constant. Modular properties then fix this t?-constant. 

Many other functional identities may be derived using the same ideas. We present here 
a (not exhaustive) list of identities, all of the general form ([3]), for the root system A 2 and B 2 
(see Appendix Bj^j. The origin of such identities stems from the analysis of elliptic solutions 
to the Witten-Dijkgraaf-Verlinde-Verlinde equation. This application will be presented in the 
next section. 

4.1 A 2 -identities 



\f^\y) + \f {l ' 2 \x + y) 



f m (x + y)[f^(x)-f^(y)] 

+f im (y)[f {2 > 1 Hx + y)-f i2 > 1 Hx)] 
+ [f [1 < 2 \x 

f^°\x)[fW(x + y)-fW( y )] 

+f {3 >°Hy)[f {1 > 2) (x + y)-f {1 > 2) (x)] 
-lf im (x + y)[f {1 > 2 \x) + f^(y)] 

+ f f {0 ' 3) (x + y) + ^E i (r) = 0. 



(7) 



+ < 



+ lf(2,l) {x + y)f {2,l) {y) 



The dependence of / on r in these formulae has been suppressed for notational convenience. 



5 Applications 

In this section we present two applications of these identities: one is in the theory of elliptic 
solutions to the Witten-Dijkgraaf-Verlinde-Verlinde (or WDVV) equations of topological quan- 
tum field theory [19], and the second, more speculative and incomplete, in the theory of elliptic 
Dunkl-type operators. 



1 The corresponding identities for G2 are available on request. 
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5.1 Frobenius manifolds and solutions to the WDVV-equations 

The basic identity ([2]) has particularly simple rational and trigonometric limits. Given a,b,c £ C 
such that a + b + c = then the rational limit gives 

111111 

-■- + -•- + -•- = 8 

a b be c a 

and the trigonometric limit gives 

cot (a) cot (b) + cot (6) cot(c) + cot(c) cot (a) = 1. (9) 

These identities are well known, and one application of them is in the construction of solutions 
to the WDVV-equation. In particular, the function 



F = \ E (a,z) 2 log(a,z)^ 



4 

satisfies the WDVV equation (see, for example, [3], H5]) 
[Fi, Fj] = 0, 

where 

d 3 F 



dz l dzidz r 



(here g is the Euclidean metric) and lZ\y are the roots of a Coxeter group W. The verification 
of this result reduces to sums over vectors IZw H 11, where il is an arbitrary 2-plane, and the 
only non-trivial configurations are the 3 irreducible 2-dimensional root systems. The WDVV 
equations are then satisfied by virtue of identities such as ([8]). Similarly, the verification that 
the trigonometric function |1CH [T6] 

F = cubic terms + h a U 3 (e i(a ' x) ) 

satisfies the WDVV equation involves use of identities of the form (J9|). It is not therefore 
surprising that the elliptic identities of the form ([2]) should play a role in elliptic solutions of the 
WDVV equation. Indeed, the origin of the differential equation in Theorem Q] is the following: 

Lemma 2. The function 

F(u, z, t) = 2 u2t ~ uz2 + r ) 
satisfies the WDVV equations if and only ifh(z,r) satisfies the partial differential equation 
h (3,0) h (l,2) _ ( /l (2,l))2 + 1^(0,3) =Q 

3 

From Theorem [TJ one solution of this equation is clearly 



2{2^if 



However, another solution is 



h(z,T)=f(2z,T)-4f(z,T) 
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and it is this solution that corresponds to a nice geometric structure, namely the dual prepoten- 
tial to the ^4i-Jacobi group orbit space [T7j. This, and the corresponding potentials that can be 
calculated for the A/v-Jacobi group orbit [T7] space suggest that one can develop the following 
functional ansatz for more general class of solutions: 

F(u, z, r) = ^u 2 r - -u(z, z) + ^2 V/((z, «), r) . (10) 

Here (u, z, r) ECx x H and it is a certain finite set of 'root vectors' in C^. Note that the 
vector field e = ^- is the unity of the multiplication, and that 

g = 2du dr — (dz, dz) , 

where ( , ) is the standard complexified Euclidean inner product on C . The precise conditions 
on the set of vectors il which are required for (|10p to satisfy the WDVV equations were derived 
in [19] . These start with the definition of a complex Euclidean V-system: 



Definition 2 ([H]). Let f) be a complex vector space with non-degenerate bilinear form (, ) 
and let il be a collection of vectors in f). A complex Euclidean V-system il satisfies the following 
conditions: 

• il is well distributed, i.e. Yl h a (a, u)(a, v) = 2/i^(u, v) for some A; 

• on any 2-dimensional plane II the set II nil is either well distributed or reducible (i.e. the 
union of two non-empty orthogonal subsystems). 

With this one can define an elliptic V-system as follows: 

Definition 3. An elliptic V-system il is a complex Euclidean V-system with the following 
additional conditions: 

• E h a (a,z) 4 = 3(z,z) 2 ; 

• the three conditions: 

hp(a,(3)(f3,a ± ) n = 0, 

/3en Q nu 

^(a,/3)(aA/3)(/3,a ± ) n = 0, 

/3en Q nu 

hp{a, f3)(a A /5) 8 (a A /?)(/?, or 1 )™ = 0; 

/3en Q nu 

• there exists a full iV-dimensional weight lattice of vectors p such that (p, a) E 7L for all 

a e il. 



With this definition of an elliptic V-system one can state the main theorem of |19j . 
Theorem 2. Let il be an elliptic V-system. If = then the function 

F(u, z, r) = ]-u 2 t - -u(z, z) + h <*f (( z > «)> r ) 



2 2 

QGll 



satisfies the WDVV equations. Ifh^^O then the modified prepotential 

F ^ F+ 3(2^F £l3(1 ' 9) 
satisfies the WDVV equations. 
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The proof of this theorem is entirely analogous to the proof of Theorem [TJ the first and third 
parts of Definition [3] ensure that that functions A = [Fj,Fj], the obstructions to the WDVV- 
equations from holding, are doubly periodic in all z- variables, and the second conditions ensure 
that these functions have no poles and hence must be a ^-constant. Modular properties then 
fix the ^-constant as above, and hence to the vanishing of the obstructions A. 

One unexpected result is the a 'pure' root system does not suffice: it is not, in general, an 
elliptic V-system. Thus generalization of classical integrable systems from rational and trigono- 
metric to elliptic based purely on Coxeter root systems may not always work. 

5.2 Dunkl-type operators 

There is a close similarity between the mathematics behind the WDVV equations and the 
commutativity of Dunkl operators. Given the role of the above functional identities in the 
construction of solutions to the associativity equations a natural question to raise is whether or 
not they have a role in Dunkl-type operators. The results presented here will be for the Coxeter 
group A2 alone. 
Let 

3(- 1 >(£) = € + X;*a(a.OSa, 

aeii 

(0 = % + Yl *^ a > 0/ (3-i '° ( z > i = 0, 1, 2, 3. 

aeil 

Here 11 = {a, (3, —(a + (3)} where a and are the simple positive roots of A%. The following 
follows from direct calculation, on using the basic identity ([2]). 

Proposition 2. Let 

V^fori) = [E( a >(0,S< 6 >fa)] + [E( 6 >(0,HWfa)]. 

Then 

F (0 '- 1) (e,r / )=0, 

(Note that the second of these can be written more succinctly as 

V [sW(0,sM(r,)]=0 

where g = dudr+2dz 2 .) These relations bears a close similarity with the zero-curvature relations 
of the elliptic WDVV equations. However, terms such as F*- ' 1 ^^, if) cannot be expressed in 
terms of other F^'- 7 ', which might have been expected via the use of the generalized Frobenius- 
Stickelberger relation (J7|). 

These operators E do have simple rational and trigonometric limits, from which one can 
recover standard results. In the rational limit /( 3 '°)(z,t) ~ z~ l and one hence obtains the 
original result [5] of Dunkl 

[S( )(£),S(%)]= 0! 

and in the trigonometric limit (q — ► 0) /^ 2,1 )(z,r) — > — i and so 

a (1) (0--^s(- 1 )(0. 
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With this one obtains Heckman's operators [13] which, following Cherednik, can be made into 
a pair of commuting operators [1] . 

It does appear, based on the above calculations for A2, that these identities cannot be used 
to construct new classes of elliptic-type Dunkl operators (thought the basic identity ([2]) does 
play a role in elliptic KZ-theory [9]). This is, however, scope for further investigation of this 
problem. 



6 Conclusions 

One recurrent theme in the theory of integrable systems is the tower of generalizations 

rational — ► trigonometric — > elliptic, 

and we have seen in this paper how the Frobenius-Stickelberger relation, via its i?-function 
form, reflects this. The WDVV-equation are themselves an example of an integrable system; 
they arise as the zero-curvature conditions of the connection [6] 

y x Y = V X Y + XX o Y, 

and the solutions described in Theorem ([2]) (at least for the Weyl groups A jy and Bjy, and 
conjecturally for all Weyl groups) sit at the right of the following tower of generalizations: 

C N /W — ► C N+1 /W — ► n/J(g) 

Coxeter group | [ Extended affine Weyl I f Jacobi group 



orbit space orbit space orbit space 

(technically the solutions are the corresponding 'almost dual' solutions [7j to the Frobenius 
manifold structures on these spaces and the solutions described in Theorem (j2j) correspond to 
the Frobenius manifolds constructed in [21 [18]). These Jacobi group orbit spaces use in their 
construction properties of Jacobi forms [U [21], functions which may be thought of as elliptic 
generalizations of ^-invariant polynomials. 

This tower of generalizations clearly need not stop at elliptic solutions. The A/v-Frobenius 
manifolds for the rational (indeed polynomial), trigonometric and elliptic cases have a Hurwitz 
space description in terms of the moduli space H g ^{ki, . . . , k{) of branched coverings of the 
Riemann sphere. An interesting question is whether or not there is an orbit space construction 
for these more general spaces: 

H ,n( n ) — H 0>N (k,N-k) — Hi iN (N) — ... — ► H g>N (N) 

III I 



C N /A N — ► C N+1 /AP — ► n/J(A N ) — f ... — ► space 



orbit 

space 
structure? 



It seems sensible to conjecture that such an orbit space exists. One would expect Seigel modular 
forms to play a role instead of the modular forms used here. Higher genus Jacobi forms certain 
have been studied, but their use has yet to percolate into the theory of integrable systems. The 
development, and applications of, the neo-classical i9-function identities studied in Section H] 
remains to be done systematically. Higher genus analogues of these identities certain exist, since 
there exists almost-dual prepotentials on these Hurwitz spaces which, by construction, satisfy 
the WDVV equations. In the genus and genus 1 cases, the prepotential is very closely related 
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to the prime form on the Riemann surface. This may be the starting point for the development 
of a functional ansatz for the higher genus cases. Central to the results presented here are 
the quasi-periodicity and modularity properties of the elliptic polylogarithm, and these were 
obtained from the analytic properties of this function; the only role the analytic properties play 
were in the development of these transformation properties. It would be attractive if one could 
obtain these directly from the geometric properties of the prime form. This approach could then 
be used in the higher genus case where the analytic properties are likely to be considerably more 
complicated. 

Appendix A 

There are, unfortunately, many different definitions and normalizations for elliptic, number- 
theoretic and other special functions. Here we list the definitions used in this paper. Let 
q = e 2mT , where r £ H. 

• -i?i-function: 

oo 

0!(z|t) = -i {e niz - e- niz ) q* JJ(1 - q n ) (l - q n e 2mz ) (l - g n e" 2 ™ 2 ) . 

n=l 

The fundamental lattice is generated by z \— > z + 1, z ^ z + t and the function satisfies 
the complex heat equation = 47rMi jT . 

• Bernoulli numbers: 

oo „ 
X „ X 



n=0 

• Eisenstein series: 

2k ^ 

E k (r) =1-—J2 °k-i(n)q n , k G 2N, 

k n=l 

where a k {n) = J2d\n d>C - 

• Dedekind ^-function: 

oo 

rj(r)=q^l[(l-q n ). 

n=l 



Polylogarithm function: 

°° z n 

U n(z) = ^2~ w , \z\<1. 



n 1 

n=l 



Note that t?i, E2 and r] are related: 
rf(r) 27ri 1 0?'(O,t) 



t](t) 24 v ' 12m #i( > r ) ' 
These have the following properties under inversion of the independent variable: 

T ~ nE n (~j =E n (r\ n> 4; 
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where in the last formula the square-root is taken to have non- negative real part. The poly- 
logarithm has the inversion property: 

(-l)*- 1 ^*- 1 ) = U N (z) + g B f^^ (logzf-i. (11) 

This may be used to analytically continue the function outside the unit disc to a multi-valued 
holomorphic function on C\{0, 1}. For a discussion of the monodromy of the polylogarithm 
function see |20|. 



Appendix B 

Here we present the identities corresponding to the -E>2-root system. There are more than one 
identity in each of these two sets: this stems from the facts that the root systems have roots of 
different lengths. 
Set B 2 (a): 

f im (x)[f^(y)-f^\x + y)] 
+ f(m {x + y )[f (2,i) {x) _ f (2,i) {x + 2y)] 

+ f^\x + 2y)[f^\x + y)-f^{y)} 



+ 



1 -f^(x)- 1 -f^(x + 2y)-2f^\y) 



f {3 '°Hy)[f i2 ' 1) (x + 2y)-f^(x)] 
+ f i3 '°Hx + y)[f^\x)-f^\x + 2y)] 
+ 2f^\x + 2y)[f^\x + y)-f^\y)] 
+ {f(^) {x)+ 2f^(x + y)-2f^(y)}=0. 

Set B 2 (b): 

/ (3,0) (a;)[/ (l,2) (a; + y) _ / (l,2) (y)] 

+ f^°\x + y)[f^(x) + f^(x + 2y)} 
+ f(m {x + 2 y)[f^(x + y)+f^(y)] 

if f^ 3 \x) + 2f^)(x + y) ) j_ 
+ 3 \ +/ (o,3) (x + 2y) + 6/ (o,3) (y) f + 3G E4iq) 

f^(y)[f {1 ' 2) (x + 2y)-f^Hx)] 
+ Af^\x + 2y)[f^(y)+f^\x + y)] 
+ f^ ) {x + y)[f^\x) + f^\x + 2y)\ 



fW(x + y)f^(x) 
+ / (2,D (x + y)/ (2,i) (x + 2y) 



f^\ x + 2y)f^ 1 \x + y) 
+ f^\x + 2y)f^(y) 



if /(0'3) (x) + 8/ (o,3) + | 

+ 3 1 +/ (0,3) (x + 2y)+8/ (0,3) (y) 



18 
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